arXiv:math/0112228v5 [math.AP] 25 Mar 2003 


WELL-POSEDNESS FOR THE LINEARIZED MOTION OF A 
COMPRESSIBLE LIQUID WITH FREE SURFACE BOUNDARY. 


Hans Lindblad 
February 8, 2008 


1. Introduction 

We consider Euler’s equations 

(1.1) p(d t + V k dk)vj = — djp, j = l,...,n in V, where di = d/dx * l , 
describing the motion of a perfect compressible fluid body in vacuum: 

(1.2) (dt + V k dk)p + pdivU = 0, divU = d k V k in V , 

where V k = 5 kl Vi = v k and we use the summation convention over repeated upper and lower indices. 
Here the velocity V = (V 1 , ...,V n ), the density p and the domain V = Uo<t<T{^} x V t , V t C M" are 
to be determined. The pressure p = p(p) is assumed to be a given strictly increasing smooth function 
of the density. The boundary dT>t moves with the velocity of the fluid particles at the boundary. The 
fluid body moves in vacuum so the pressure vanishes in the exterior and hence on the boundary. We 
therefore also require the boundary conditions on dfD = Uo<t<T{^} x dT> t : 

(1.3) {d t + V k d k )\ av sT(&D), 

(1.4) p = 0, on dV. 

Constant pressure on the boundary leads to energy conservation and it is needed for the linearized 
equations to be well posed. Since the pressure is assumed to be a strictly increasing function of the 
density we can alternatively think of the density as a function of the pressure and for physical reasons 
this function has to be non negative. Therefore the density has to be a non negative constant ~p 0 on the 
boundary and we will in fact assume that p 0 >0, which is the case of liquid. We hence assume that 

(1.5) p(p 0 ) = 0 and p'(p) > 0, for p > p 0 , where p 0 > 0 

From a physical point of view one can alternatively think of the pressure as a small positive constant 
on the boundary. By thinking of the density as function of the pressure the incompressible case can be 
thought of as the special case of constant density function. 

The motion of the surface of the ocean is described by the above model. Free boundary problems 
for compressible fluids are also of fundamental importance in astrophysics since they describe stars. 
The model also describes the case of one fluid surrounded by and moving inside another fluid. For large 
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massive bodies like stars gravity helps holding it together and for smaller bodies like water drops surface 
tension helps holding it together. Here we neglect the influence of gravity which will just contribute 
with a lower order term and we neglect surface tension which has a regularizing effect. 

Given a bounded domain Vq C R n , that is homeomorphic to the unit ball, and initial data Vo and 
Po, we want to find a set T> C [0, T] x R n , a vector field V and a function p, solving (1.1)-(1.4) and 
satisfying the initial conditions 

(1.6) {x; (0,x) G T>} = T> 0 , 

(1.7) V = V 0 , p = po on {0} x Vo- 

In order for the initial-boundary value problem (1.1)-(1.7) to be solvable initial data (1.7) has to satisfy 
certain compatibility conditions at the boundary. By (1.2), (1.4) also implies that &rvV\ dv = 0. We 
must therefore have Po\ aT , 0 = 7 >o and div Vq| gT> ^ = 0. Furthermore, taking the divergence of (1.1) gives 

an equation for (dt + V k dk) divV in terms of only space derivatives of V and p, which leads to further 
compatibility conditions. In general we say that initial data satisfy the compatibility condition of order 
m if there is a formal power series solution in t, of (1.1)-(1.7) (p,V), satisfying 

(1-8) (dt + V k dky(p-p 0 )\ {0}xaVo = 0, j = 0, ..,m — 1 

Let J\f be the exterior unit normal to the free surface dT> t . Christodoulou[C2] conjectured the initial 
value problem (1.1)-(1.8), is well posed in Sobolev spaces under the assumption 

(1.9) V/vp < —c 0 < 0, on dV , where \7j^=AT l d x i. 

Condition (1.9) is a natural physical condition. It says that the pressure and hence the density is larger 
in the interior than at the boundary. Since we have assumed that the pressure vanishes or is close to 
zero at the boundary this is therefore related to the fact that the pressure of a fluid has to be positive. 

In general it is possible to prove local existence for analytic data for the free interface between two 
fluids. However, this type of problem might be subject to instability in Sobolev norms, in particular 
Rayleigh-Taylor instability, which occurs when a heavier fluid is on top of a lighter fluid. Condition (1.9) 
prevents Rayleigh-Taylor instability from occurring. Indeed, if this condition is violated Rayleigh-Taylor 
instability occurs in a linearized analysis. 

In the irrotational incompressible case the physical condition (1.9) always hold, see [W1,2,CL], and 
[Wl,2] proved local existence in Sobolev spaces in that case. [Wl,2] studied the classical water wave 
problem describing the motion of the surface of the ocean and showed that the water wave is not 
unstable when it turns over. Ebin[El] showed that the general incompressible problem is ill posed in 
Sobolev spaces when the pressure is negative in the interior and the physical condition is not satisfied. 
Ebin[E2] also announced a local existence result for the incompressible problem with surface tension on 
the boundary which has a regularizing effect so (1.9) is not needed then. 

In [CL], together with Christodoulou, we proved a priori bounds in Sobolev spaces in the general 
incompressible case of non vanishing curl, assuming the physical condition (1.9) for the pressure. We 
also showed that the Sobolev norms remain bounded as long as the physical condition hold and the 
second fundamental form of the free surface and the first order derivatives of the velocity are bounded. 
Usually, existence follows from similar bounds for some iteration scheme, but the bounds in [CL] used 
all the symmetries of the equation and so only hold for modifications that preserve all the symmetries. 
In [LI] we showed existence for the linearized equations and in [L2] we proved local existence for the 
nonlinear incompressible problem with non vanishing curl, assuming that (1.9) holds initially. 
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For the corresponding compressible free boundary problem with non-vanishing density on the bound¬ 
ary, there are however in general no previous existence or well-posedness results. Relativistic versions 
of these problems have been studied in [C1,DN,F,FN,R] but solved only in special cases. The methods 
used for the irrotational incompressible case use that the components of the velocity are harmonic to 
reduce the equations to equations on the boundary and this does not work in the compressible case 
since the divergence is non vanishing and the pressure satisfies a wave equation in the interior. To be 
able to deal with the compressible case one therefore needs to use interior estimates as in [CL,LI]. Let 
us also point out that in nature one expects fluids to be compressible, e.g. water satisfies (1.5), see [CF], 
For the general relativistic equations there is no special case corresponding to the incompressible case. 

Here we show existence for the linearized equations and estimates for these in Sobolev spaces in the 
general compressible case (l.l)-(l.T), assuming that (1.8) and (1.9) hold. This can be considered as a 
linearized stability result, showing that small perturbations of initial conditions in Sobolev spaces leads 
to small perturbations for finite times. Furthermore, in a forthcoming paper [L3] we use existence and 
estimates for the inverse of the linearized operator to prove existence for the nonlinear problem using 
the Nash-Moser technique also in the compressible case. The existence proof here uses the orthogonal 
decomposition of a vector field into a divergence free part and a gradient of a function that vanishes 
on the boundary. For the divergence free part we get an equation of the type studied in [LI] and for 
the divergence we get a wave equation on a bounded domain with Dirichlet boundary conditions. The 
interaction terms between these equations are lower order so if we put up an iteration, the equations 
decouple for the the new iterate and the previous iterates only enter in the lower order terms. 

Existence of solutions for the wave equation on a bounded domain is well known. However dealing 
with the divergence free part of the equation requires the techniques developed in [LI]. Here we use a 
generalization of the existence theorem in [LI] to the case when the divergence of the solution we linearize 
around is non vanishing. In [LI] we showed that the linearized incompressible Euler’s equations becomes 
an evolution equation for what we called the normal operator. The normal operator is unbounded and 
not elliptic in the case of non vanishing curl. It is however positive assuming the physical condition (1.9) 
and this leads to existence. Up to lower order terms, the projection of the linearized compressible Euler’s 
equations onto divergence free vector fields becomes the linearized incompressible Euler’s equations. 

As pointed out above, the positivity of the pressure (1.9) leads to the positivity of the normal 
operator, introduced in [LI]. It appears that this condition is needed for the well-posedness also in the 
compressible case since the divergence free part essentially decouples from the divergence. In fact, the 
compressible case was the main motivation for formulating (1.9) since in that case it is clear that the 
pressure has to be positive and in nature one expects fluids to be slightly compressible. 

In order to formulate the linearized equations one has to parametrize the boundary. Let us therefore 
express Euler’s equations in the Lagrangian coordinates given by following the flow lines of the velocity 
vector field of the fluid particles. In these coordinates the boundary becomes fixed. Given a domain Dq 
in R", that is diffeomorphic to the unit ball Q, we can by a theorem in [DM] find a diffeomorphism /o : 
D —> Vo with prescribed volume form det ( dfo/dy ) up to a constant factor. Let V and vGC(V) satisfy 
(1.3). The Lagrangian coordinates y are given by solving for the Eulerian coordinates x = x(t, y) = ft(y) 


( 1 . 10 ) 


dx(t, y) 
dt 


= V ( t , x(t, y)), x(0, y) = f 0 (y), y € Q. 


Then f t : D —> V t is a diffeomorphism and the boundary becomes fixed in the new y coordinates. Let 


( 1 . 11 ) 




T rk d d 

+ v L -— = — 

x=const dx k dt 


and 


d„ = 


d dy a d 


y=const 


dx l dx l dy a ’ 


be the material derivative and partial differential operators expressed in the Lagrangian coordinates. 
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In these coordinates Euler’s equation (1.1), the continuity equation (1.2) and the boundary condition 
(1.4) become: 

(1.12) D 2 x l + dih = 0, D t p + p divF = 0, p\ d n = Poi 


where the enthalpy h = h(p) = jf p'(p)p 1 dp is a strictly increasing function of p 


and x, V = D t x 


and p are functions of (t, y) € [0, T] X Q. Furthermore, p can be determined from x: 

(1.13) p = khi~ 1 , where k = det (dx/dy) and k = pn | , 

since D t k = /-cdivH. In (1.10) there is a choice of mapping f 0 and domain fh By [DM] one can find a 
diffeomorphism with prescribed volume form up to a constant between any two diffeomorphic sets. We 
therefore choose D to be the unit ball and det (dfo/dy) so det ( dfo/dy)po = k is any given fixed function 
k(y) that we take to be constant. Making this choice, initial data for p is part of the initial data for x. 

The free boundary problem for Euler’s equations (1.1)-(1.7), hence become an equation for x(t,y ): 


-l 


(1.14) D/x 1 + &i h = 0, p = kn , 


m 


[0, T] xH, n = 1, where d- L = 


8y a d 


I aO ’ 1 Q x i Qya ’ 

h=h(p) is a strictly increasing function of pand p=p(n) is a function of k = det (dx/dy). Initial data are 


(1.15) 


x 


t =0 


= /o, D t x \t=o = ^0 


In order for (1.14) to be solvable, initial data has to satisfy the constraints; det (dfo/dy) | = 1, 

divVo| 0 Q = 0 and taking the divergence of (1.14) gives an equation for D t divV in terms of space 
derivative of x and V = D t x which leads to further conditions. Since (1.14) gives D/x in terms of space 
derivatives of x we can obtain a formal power series solution in time t, x, to the first two equations in 

(1.14) satisfying the initial conditions (1.15). The compatibility condition of order m is the requirement 
that the formal power series solution up to terms of order m satisfy the boundary condition in (1.14): 

(1.16) D J t (det (dx/dy) - l)| Qxan = 0, j = 0, ...,m-l 


Let us now derive the linearized equations. (1.14) can be thought of as an equation <h(x) = 0, where 
is a functional of x(t,y) given by $(x)* = D^x 1 + d t h. for 1 < i < n, where h is a given function of 
k = det (dx/dy) and d % are the differential operators in (1.14) with coefficients depending on derivatives 
of x as well, and 4>(x) n _|_i = (k — l)| 0fr We assume that x(t,y) is a given smooth solution of (1.14), 
i.e. <I > (.t) = 0. Let x(t,y,r) be a smooth function also of a parameter r, such that x\ r _ Q = x and set 
5x = dx/dr | Q . Then the linearized equations are the requirement on 6x, that x satisfies the equations 

(1.14) up to terms bounded by r 2 as r —» 0, i.e. 4>'(x)((5x) = d$(x)/dr | = 0. If we replace x in 

(1.14) by x and apply 6 = <9/<9r|^ =0 we hence obtain the linearized equations: 

(1.17) D/5x l + (didkh)5x k — di((dkh)Sx k — Sh ) = 0, 5h = —h'(p)p divdx, div<5x| an = 0. 

Here we used that [<5,<9,;] = —(di5x k )dk and Sp = ~pk~ 1 5k = —pdivdx, see section 2 and [LI]. The 
initial data for the linearized equations are 

(1.18) Ht=o = S f°> Dt Sx \t=o = SV °‘ 

The initial data are as before subject to constraints. Let 5x be the formal power series solution in time 
t to (1.17)-( 1.18). The compatibility condition of order m is 

(1.19) D\ divhx = 0, j = 0,..., m — 1. 

The main difference between (1.17) and (1.14) is the higher order term di((dkh)5x k ) , since the term 
didh , depending on divdx, in (1.17) corresponds to the term dih, depending on det (dx/dy) in (1.14). If 
we take x above to be a family of solutions of (1.14) depending on the parameter r, then our estimates 
below show that a small change of initial conditions only give rise to a small change of the solution in 
Sobolev spaces. Our main result is the following linearized stability result: 
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Theorem 1.1. Let Q be the unit ball in R n and suppose that x is a smooth solution of (1.14) satisfying 
(1.9) for 0 < t < T. Suppose that (8fo,8Vo) are smooth satisfying the compatibility conditions of all 
orders m, i.e. (1.19) holds for all m. Then the linearized equations (1.17) have a smooth solution 
Sx for 0 < t < T satisfying the initial conditions (1.18). Let Af be the exterior unit normal to dV t 
parametrized by x(t,y) and let 8xjy = Af ■ 8x be the normal component of 8x. Set 

(1.20) E r (t) = || D t 8x(t, •)||if'-(n) + II 8x(t, ■)l|/f r (n) + I! di v8x(t, -)llff'-(n) + \\$ x N(t, ‘)llH r ( 9 n) 

where H r {Ll) and H r (dLl ) are the Sobolev spaces in fl respectively on dLl. Then there are constants C r 
depending only on x, r and T such that 

(1.21) E r (t) < C r E r ( 0), for 0 < t < T, r > 0. 

Furthermore, let N r (Q) be the completion of C°°(Q) in the norm ||&c||.H>(fi) + || div&c(i, •)l|// r (n) + 
||<fo/v||.H>(^)- Then if the constraints in (1.19) hold for all orders m and 

( 1 . 22 ) (8f 0 ,8V 0 )€N r (n)xH r (n) 

it follows that (1.17)-(1.18) has a solution 

(1.23) 0 5x,D t 8x ) G C([0,T],N r (Q) x H r (Q)). 


As we have argued, any smooth solution of (l.l)-(l.T) with Vo diffeomorphic to the unit ball can be 
reduced to a smooth solution of (1.14) where Ll is the unit ball. That there are initial data (1.18) such 
that (1.19) hold for all m follows by taking S fo and 8Vq compactly supported in the interior of fl. The 
term ||&E/v'l|.H’ T ’(0n) is equivalent to the variation of the second fundamental form 6 = dAf of the free 
boundary dV t measured in H r ~ 2 (dLl). For a general component we can only say that 8x£H r ~ 1 / 2 (dLl). 
The energy estimate (1.21) also hold in the incompressible case when div<5x = 0, see [LI], and in [CL] 
we obtained similar bounds for ||u||#r(n) + ||0 ||h»— 2 (am in the nonlinear incompressible case. 

Let us now outline the main ideas in the proof. We will rewrite the linearized equations (1.17) in a 
geometrically invariant way and use this to obtain energy bounds and existence. We have defined our 
vector fields as functions of the Lagrangian coordinates (t, y) G [0, T\ x LI but we can also think of them 
as functions of the Eulerian coordinates (f, x) G V, and we will make this identification without saying 
that we compose with the inverse of the change of coordinate y —> x(t,y). The time derivative has 
a simple expression in the Lagrangian coordinates but the space derivatives have a simpler expression 
in the Eulerian coordinates, see (1.11). For the most part we will think of our functions and vector 
fields in the Lagrangian frame but we use the inner product coming from the Eulerian frame, i.e. in 
the Lagrangian frame we use the pull-back metric of the Euclidean inner product: 

(1.24) X • Z = SijX'Z* = g ab X a Z b , where X a = X'-^~, g ab = Sij ——. 

Here X 1 refers to the components of the vector X in the Eulerian frame, X a refers to the components 
in the Lagrangian frame, g a b is the metric in the Lagrangian frame and 8j 3 is the Euclidean metric. The 
letters a, b, c, d, e, /, g will refer to indices in the Lagrangian frame whereas i, j, k, l, m, n will refer to the 
Eulerian frame. The norms and most of the operators we consider have an invariant interpretation so 
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it does not matter in which frame they are expressed. In the introduction we express the vector fields 
in the Eulerian frame but later we express them in the Lagrangian frame. The L 2 inner product is 

(1.25) {X,Z)= [ X-Zdx= [ X • Z ndy = [ X ■ Z p” 1 k(y) dy , 

Jv t JQ. J f! 

where k = det (dx/dy) and k = Kp | . 

Let us first point out that the boundary condition p\ gn = Po leads to that the energy is conserved 
for a solution of Euler’s equations (1.12). Let Q{p) = f-f 2q(p)p~ 2 dp, where q(p) = p(p) — p(j> 0 ). Since 
D t (pK) = 0 and pD t Vi = —dip(p) = —diq(p) it follows from the (1.25) and the divergence theorem that 

(1-26) li f d y l 2 + Q{p))pdx = f D t (\V\ 2 + Q(p))pdx = [( - 2V l diq(p) + 2q(p)p~ 1 D t p) dx 

at JT> t JT> t JV t 

= 2 [ (divE q(p) + q(p)p~ 1 D t p) dx - 2 f Vxq{p) dS = 0, 

Jv t JdD t 

where Vj\f= A/)E* is the normal component and we also used (1.12). 

We will obtain similar energy estimates for the linearized equations (1.17) for energies containing 
an additional boundary term. We will first rewrite the linearized equations in a geometrically invariant 
way. The last term in the first equation in (1.17) is a positive symmetric operator in the energy inner 
product on vector fields satisfying the boundary condition divX\ dv ^ = 0: 

(1.27) CX = -V(h'(p)(p divX + ( d k p)X k )) = -V(h’{p) div(pX)), where V* = S l] d J 

i.e. (X,pCX) > 0 and (Z,pCX) = (pCZ,X) if &\vX\ dT) * = di vZ\ dT> ^ = 0. In fact, if divX | dT) ^ = 0 
then h'(p) div(pX)\ gT> ^ = X k df-h \ gv = XnP~ 1 XnP\ g' Dt an< ^ integrating by parts we get 

(1.28) ( Z,pCX)=[ div(pZ) div(pX)h\p)dx+ ( Z N X N (-X N p)dS 

Jv t JdV t 

which proves the symmetry and the positivity follows since h'(p) > c\ > 0 and VnP < — Co < 0. 

We will also replace the time derivative by a time derivative that preserves the boundary condition. 
Let 

(1.29) C D ,X‘ = D t X‘ ~ (9i-V)X‘ = —D^X*), 

be the space time Lie derivative with respect to D t = (1, V) restricted to the space components or 
equivalently, the time derivative of the vector field expressed in the Lagrangian frame. Let 

(1.30) C D ,X< = C D ,X‘+ divCX* = k -i^D,( k |L.X‘) 

be the modified Lie derivative that preserves the boundary condition, divX| af2 = 0. In fact 

div£ Dt X = D t divX, where D t = D t + div V. 

6 


(1.31) 



The linearized equations (1.17) can now be written as an evolution equation for the operator C: 


(1.32) A + CX = B(X,X), divX| t=o = 0, 

where X = Sx, X = Co t X, X = C? D X and B is a linear form with coefficients depending on x and p. 
Associated with (1.32) is the energy 

(1.33) E(t) = (X,pX) + (X,p(C + I)X), 

and we prove that E'<CE which gives the bound (1.21) for r= 0. The boundary term in (1.20) comes 
from (1.28). To obtain estimates for higher order derivatives one can apply modified Lie derivatives with 
respect to tangential vector fields as in [LI]. This does however not prove existence for (1.32) which is 
non-standard since C is non-elliptic, time dependent and the boundary condition is non-trivial. 

We use the orthogonal projection onto divergence free vector fields in the inner product (1-25) to 
obtain an equation for the divergence and an equation for the divergence free part. The equations 
decouples to highest order, and existence and estimates for the system follows from existence and 

estimates for each equation with an inhomogeneous term. The orthogonal projection is 

(1.34) PX = X — Vq, where A q = div X, q =0. 

We will obtain a system of equations for Xq = PX and X\ = (/ — P)X by projecting the linearized 

equation (1.32) onto divergence free vector fields respectively the orthogonal complement. Taking the 
divergence of (1.32) gives a wave equation for div A with Dirichlet boundary condition: 

(1.35) Df divW - A(h'{p)pdivX 1 ) = A(X k d k h) +di .vB(X,X), divA^ = 0 

for which existence is known if h'(p)p and the metric g a bj hidden in A = )T)” =1 df = d a ng ab db , are 
bounded from above and below and the right hand side is thought of as a known function, see section 
6. Ai is then determined from divAi by solving the Dirichlet problem: 

(1.36) Ai = V< 7 i, Ac/i = divAi, q\ =0. 

To obtain an equation for the divergence free part Aq we project (1.32) onto divergence free vector 
fields. It follows from (1.27) that 

(1.37) AX = PCX = P(-V(X k d k h)), 

since div A | = 0 and the projection of a gradient of a function that vanishes on the boundary vanishes 
The operator A is a positive symmetric operator on divergence free vector fields, if condition (1.9) holds: 

(1.38) (X,AZ) = - [ X i d i {Z k d k h)dx = [ X N Z N (-X N h)dS, if divA = divZ = 0. 

Jv t JerD t 

Furthermore, we note that the commutator of time derivatives with the projection is lower order 

[C Dt ,P]X 1 = -P{(£ Dt 5 ij )6 jk (I - P)X k ) 
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which follows since Co t preserves the divergence free condition and the projection of 8 l ^djD t q vanishes 
if q\ gn = 0 and hence D t q\ ;)Q = 0. Hence PX\ = PB- 2 (X i,Xl) can be determined in terms of X\ and 
X\. Projection of (1-32) therefore gives an evolution equation for the operator A for the divergence free 
part: 

(1.40) X 0 + AX 0 = -PB 2 (X 1 ,X 1 )-AX 1 +PB(X,X) 

Existence for (1.40) with the right hand side thought of as a known function is a generalization of the 
existence proof in [LI]. For the divergence part we have an equation which is equivalent to (1.35)-(1.36): 

(1.41) X 1 -V(h\p)pdivX 1 )-PB 2 (X 1 ,X 1 ) = (. I-P)V((d k h)X k )+(I-P)B(X,X ), divAT 1 | 0O = 0 

We will show existence for the system (1.40)-(1.41) for (Xo,Xi) and from that we obtain a solution 
X = X 0 + Xi to (1.32), since these equations are exactly the projection of (1.32) onto the divergence 
free vector fields respectively the orthogonal component. The system (1.40)-(1.41) can be solved by 
iteration. If X is an iterate, then from (1.35)-(1.36) and (1-40) we get X\ and Xq and a new iterate is 
Xq+X\. There is no loss of regularity in this procedure since divA has the same space regularity as X. 

However, in order for it to be possible to solve (1.35) the initial conditions and the equation must be 
compatible with a formal power series solution satisfying the boundary conditions, and we must make 
sure that this is true at each step of the iteration. (1.32) gives X in terms of only space derivatives of 
X and X and we hence obtain a formal power series solution in time, the first two terms coming from 
the initial conditions. Our assumption is that this formal power series solution satisfies the boundary 
condition. From this formal power series one can construct an approximate solution X satisfying the 
initial conditions, the equation to all orders as t —* 0, and the boundary condition. We can then take the 
approximate solution as our first iterate or equivalently subtract off the approximate solution from A, 
which produces an inhomogeneous term vanishing to all orders as t —> 0 and vanishing initial conditions. 

Let, us now conclude the introduction by giving the main estimates we use. Since the time derivative 
preserves the boundary condition it is natural to use norms which also contain time derivatives up to 
full order in the proof, and the estimate in Theorem 1.1 afterwards follow from these. Let 

(1.42) \\X(t)\\ H r = \\X(t, 011^(0), ||A(t)|| r = 2 \\ D tX(t)\\Hs, (Xo(t)) r = ||A 0JV (i, Olltf-on) 

s-\-k<r 

where Aon = Xq ■ N is the normal component. For the divergence free equation: 

(1.43) A 0 + HA 0 = F 0 , PF 0 = F 0 
we have the estimate 

(1.44) j|A 0 (t)|| r + ||A 0 (t)|| r + (X 0 (t)) r < c(||Xo(0)|| r + ||A o (0)|| r + {X o (0))r + J 11^0Hr dr). 

This is a generalization of the estimate for the incompressible case in [LI]. For r = 0, one uses the 
symmetry and positivity (1.38) of A to prove that E = (X 0l X 0 ) + (X 0 , (A + I)X 0 ) satisfies E' < CE. 
Note that the boundary term comes from using (1.38). For r > 0, it follows from commuting modified 
Lie derivatives with respect to tangential vector fields through the equation to obtain similar equations 
and estimates for these, together with better estimates for the curl since the curl of A vanishes. For 

(1.45) A\ - V(h , (p)pdWX 1 ) - PB 2 (X 1 ,X 1 ) = Tj , (I - P)F 1 = F u divA,^ = 0, 
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we have the estimate 


(1.46) ||X 1 (i)|| r+1 < c (\\ X 1 ( 0)\\ r+1 + J ||Fi||r-idr). 

The last estimate follows from estimating the wave equation (1.35) with the right hand side replaced 
by divT\ and inverting the Laplacian (1.36). In (1.46) we do not need space derivatives up to highest 
order of F\ , since one obtains space derivatives from time derivatives through inverting the Laplacian 
in the wave equation. Using that the right hand side of (1.41) is a gradient (1.46) also holds for r = 0. 

With F 0 equal to the right hand side of (1.40) and F\ equal to the right hand side of (1.41) the 
norms in the integrals in (1.44) and (1.46) can be estimated by the sum of the norms in the left of (1-44) 
and (1.46) and this gives a priori bounds as well as uniform estimates for iterates, if we let the Fq 
and Fi be obtained from the previous iterate and solve (1-43) and (1.45) for the new iterate. Once we 
obtained the solution to the system (1.40)-(1.41), X = Xq + X\ is the solution to (1.32). The norm in 
(1.20) is bounded by the sum of the norms in the left of (1-44) and (1.46). Using (1.32) one can bound 
time derivatives in terms of space derivatives and using that the projection is continuous in the norms 
(1.42) it follows that the norms in the left of (1.44) and (1.46) can be bounded by (1.20). 


2. Lagrangian coordinates and the linearized equation. 


Let us introduce Lagrangian coordinates in which the boundary becomes fixed. Let 14 be a domain 
in R” and let /o : 14 —> Vq be a diffeomorphism. We assumed that Vq is diffeomorphic to the unit 
ball and that v(t,x), p(t,x), (t,x) £ V are given satisfying the boundary conditions (1.3)-(1.4). The 
Lagrangian coordinates y are given by solving for the Eulerian coordinates x = x(t,y) = ft(y) in 


( 2 . 1 ) 


dx/dt = V(t,x(t,y)), x(0,y) = f 0 (y), y £ 14 


Then f t : 14 —> V t is a diffeomorphism, and the boundary becomes fixed in the new y coordinates. Let 
us introduce the notation 


( 2 . 2 ) 


D ' = i 


y=constant 


d_ 

di 


+ V k ——^ • 

x=constant dx k 


for the material derivative and 
(2.3) 


di = 


d dy a d 


dx i dx i dy a ' 


for the partial derivatives. 

In these coordinates Euler’s equation (1.1) become 


(2.4) pDjxi + dip = 0, 
and the continuity equation (1.2) become 

(2.5) D t p + p divU = 0, 


(t,y) £ [0, T\ x 14 


(t,y) £ [0, T] x 14 


Here the pressure p = p(p) is assumed to be a given smooth strictly increasing function of the density 
p. Let p 0 be defined by p(~p 0 ) = 0. Let h. the enthalpy, be defined by 


( 2 . 6 ) 


Hp) = [_ p'{p)p 1 dp. 
dp 0 



Then (2.4) becomes 

(2.7) 


D^Xi + dih = 0, (t, y ) G [0, T] x D. 


The density p satisfies (2.5) but since 

(2.8) re = det (dx/dy), 
satisfies 

(2.9) D t K — K divH = 0 

it follows that p = pqKq/k, where po and kq are the initial values. By a theorem of [DM] one can 
arbitrarily prescribe the volume form reo up to a constant so we take reo = k/po, where A: is a constant, and 
Q to be the unit ball, by composing with a diffeomorphism, since we assumed that T >o is diffeomorphic 
to a unit ball. Hence p is determined from x: 


(2.10) p = /c/re 

By choosing the constant k appropriately the boundary condition (1.4) can hence be expressed 


( 2 . 11 ) 


lan 


= 1 


Since h is a function of p which in turn by our choice (2.10) is a function of re = det (dx/dy) we can 
think of h as a function of re. (2.7) is then an equation involving the coordinate x only and initial data 
for p is included in the choice of initial mapping /q. Initial data for (2.7) are 

(2-12) ar | t=0 = /o, D t x\ t=0 = V 0 

In order for (2.7) to have a smooth solution satisfying (2.11), initial data has to satisfy the constraints 
det (dfo/dy)\ gc = 1 and divVo| af2 = 0, by (2.9). Taking the divergence of (2.7) gives 

(2.13) D t divV + (d i V k )d k V i + Ah = 0 

which leads to further conditions. Since (2.7) gives D/x in terms of space derivatives of x we can 
obtain a formal power series solution in time t. x, to (2.7) satisfying the initial conditions (2.12). The 
compatibility condition of order m is the requirement that the formal power series solution up to terms 
of order m satisfy the boundary condition in (2.11): 

(2.14) D J t ( det (dx/dy) - l)| Qxan = 0, j = 0,...,m-l 


At this point we also remark that we get a wave equation for h. Since h is a strictly increasing 
function of p we can think of p = p(h) as function of h. Hence with e(h) = In p(h) (2.5) instead become 


(2.15) D t e(h) + divH = 0. 

and this together with (2.13) gives a wave equation for h with Dirichlet boundary conditions: 


D 2 t e(h) - Ah - (d i V k )d k V i = 0, 
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= 0 . 

on 


(2.16) 


h 



Here 


(2.17) A h = Y^dfh = n~ l d a (i%g ab d b h), where g a b = hj^^, 

is the metric in the Lagrangian coordinates and g ab is its inverse. Here d a = d/dy a and we use the 
convention that differentiation with respect to the Eulerian coordinates is denotes by letters i , j, k , l,m,n 
and with respect to the Lagrangian coordinates is denoted by a,b,c,d,e, f . In order for (2.16) to be 
solvable we must have that 

n 

(2.18) 0<e / + l/e / <ci ^ {\g ab \ + \g a b\) < n cf , \dx/dy\ 2 + \dy/dx\ 2 < c\ 

a,6=1 

for some constant 0 < c\ < oo. The first condition is related to that the pressure is assumed to be a 
strictly increasing smooth function of the density. The second and third condition are equivalent and 
says that the coordinate mapping is a diffeomorphism. Furthermore, it is well-known that one needs 
compatibility conditions to solve (2.16). 

Let us now derive the linearized equations. The calculations that follows below are similar to those 
in [LI] since the equation (2.7) mathematically is the same as the equation for the incompressible case 
with the enthalpy h replaced by the pressure p. We therefore refer the reader to [LI] for more details. 

We now assume that we have a smooth solution x = x(t,y) of (2.7) satisfying (1.9) for 0 < t < T 
and we will derive the linearized equations at this solution. Assume that x = x(t, y, r ) is a smooth 
function also of the extra parameter r such that x| 0 = x and set 5x = dx/dr | r _ 0 - Then the linearized 
equations are the requirements on 5x that x satisfies (2.7) and (2.10)-(2.11) up to terms bounded by r 2 
as r—y 0. Let 5 be a variation in the Lagrangian coordinates, i.e. a derivative 

(2-19) Sf = df/dr\ r=0 . 

with respect to the parameter r when t and y are fixed. Then [h, D t \ = 0, 

(2.20) [<S, di] = -(diSx k )d k , 
so [h— 5x k dk, di] =0. Applying S — 5x k d k to (2.7) gives: 

(2.21) D 2 5xi — ( dkD 2 Xi)5x k — di(6x k dkh — 8h ) = 0 

Since h = h(p) where p = k/n and 

(2.22) 5k = ftdivfe. 
it follows that 

(2.23) Sh = — h'(p)pdiv5x. 

The variation of the boundary condition (2.11) become 

(2.24) divhx| 0fi = 0. 
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The initial data for (2.21) with Sh given by (2.23) are 

(2.25) Sx = Sf 0 , D t Sx = SV 0 

In order for it to be possible to have a smooth solution of (2.21) and (2.23)-(2.24) initial data (2.25) must 
satisfy certain compatibility conditions. The initial data are subject to the constraints divd/o| sn = 0 
and divd'Vo = (diV k )dk5x k . 

Taking the divergence of (2.21) using (2.23) and (2.24) gives a wave equation for div<5x with Dirichlet 
boundary conditions: 

(2.26) D 2 div<Sx - Sx'diDt divV - A(5x k d k h - Sh ) + 2 (d i V k )d k (5V i - Sx^V*) = 0, 

which gives further conditions. Since (2.21) gives D 2 5x in terms of only space derivatives of 5x, this 
gives a formal power series solution in time t, which we call Sx. The compatibility condition of order 
m is the requirement that the formal power series solution satisfies the boundary condition (2.24): 

(2.27) D J t div&c = 0, j = 0, ...,m — 1. 

The basic assumption in solving the system (2.21)-(2.25) is that one should assume that div<fx has the 
same space regularity as Sx. 

Let us now express also the vector field in the Lagrangian frame. Let 

dy a 

(2.28) r = ^ 

Then, 

(2.29) D t 5x i = D t (lT b <9x7<9y 6 ) = (D t W b )dx i /dy b + WhdV^dy 1 = (D t W b )dx i / dy b + Sx k d k V i 
and multiplying with the inverse dy a /dx l gives 

(2.30) D t W a = ^-C Dt Sx\ and D t W a = ^-Ld.Sx 1 . 

ax % ox 1 

where the Lie derivative and modified Lie derivative are given by (1.29)-(1.30) and 

(2.31) D t W a = D t W a + (divV)W a = k” 1 D t {nW a ). 

Since the divergence is invariant 

(2.32) divfe = divlT = 
it therefore follows that 

(2.33) div D t W = D t divW. 

Differentiating (2.30) once more gives 

D 2 Sx i - ( d k D t V i )5x k = {D 2 W b )dx i /dy b + 2{D t W b )dV i / dy b 
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(2.34) 



It follows that 


(2.35) 


W‘ {D ‘ Sx ' - {aMV ' )Sxi) = WW^ wt + wvt) WW 8 ‘ v ’ 

= g ab D 2 W b + (A g ab - u) ab )D t W b 


where g ab is given by (2.17) and 

c)t^ r)T^ r) r r^ 

(2.36) D t g a b = {diVj + djVi) , o; ab = — (diVj - djVi). 

With d a = d/dy a the linearized equation (2.21) and (2.23) become 

(2.37) g ab D 2 t W b - d a ((d c h)W c - Sh ) = ~(D t g ac - tu ac )D t W c , 5h = -// divW 

Let Dt be as in (2.31), i.e. D t = ( D t + a), where a = In k and & = D t a = divIL Then 

(2.38) D 2 = D 2 — 2&D t + a 2 — a, D t = D t — &. a = D 2 a. 

Hence, with W = D t W and W = D 2 W , we can write (2.37) as LW = 0, where 

(2.39) LW = W a - g ab d b ((d c h)W c - 5h) - B a {W,W), 5h = -p' divW 
where 

(2.40) B a (W,W) = -g ab (g bc - to bc ){W c - &W C ) + 2 &W a + (a - a 2 )W a . 

3. The compatibility conditions, statement of the 

THEOREM AND THE LOWEST ORDER ENERGY ESTIMATE. 

We now consider the linearized operator 

(3.1) LW = W+ CW - B(W,W) 

where W = DtW, W = D 2 W , D t = D t + (divF), B is the bounded operator given by (2.40) and 

(3.2) CW a = -g ab d b ({d c h)W c +p'dwW). 

We want to show existence and estimates for the linearized equations with an inhomogeneous term F, 

(3.3) LW = F, 
with initial data 


n =0 =^o, h=o=^ 

divTy L = °- 
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(3.4) 

and boundary data 

(3.5) 



The reason for the inhomogeneous term F is that one can reduce to the case of vanishing initial 
data and an inhomogeneous term F that vanishes to all orders as t —> 0 and it is easier to first prove 
existence for this case. Differentiating (3.3) with respect to time we get 

(3.6) D k+2 W = B k (W ,., D k+1 W,dW ,..., dD k W, d 2 W, ...,d 2 D k W) + D k F, 
for some function B k . Let us therefore define functions of space only by 

(3.7) W k+2 = B k (W 0 ,...,W k+1 ,dW 0 ,... 1 dW k ,d 2 W 0 ,...,d 2 W k )\ t=0 + D^F\ t=0 , k> 0 
In view of (3.5) it follows that 0 = D k divW\ dn = div D k W\ dn so we must have 

(3.8) divIT fc | 0 ^ = 0, k = 0,...,m 

(3.8) is called the m th order compatibility condition and in order for it to be possible for (3.3)-(3.5) to 
have a smooth solution these have to hold for all orders m. 

We now define the approximate power series solution by 
(3-9) W(t,y) = ^2x(t/£k)Wk(y)t k /kl. 

Here x is smooth x( s ) = 1 for |s| < 1/2 and x( s ) = 0 f° r l s l ^ 1- The sequence e k > 0 can be chosen 
so that the series converges in C m ([ 0, T], H m ) for any m if take (HWfcll^fe + \)e k < 1/2. It follows that 

(3.10) divir| af2 = 0. 

(3.9) multiplied with n(t, y ) is a power series expansion of kW and it hence follows that (3.9) satisfies 
(3.3)-(3.5) to all orders as t —» 0: 

(3.11) D k t (LW-F) | t=() = 0, k = 0,... 

It follows that we can reduce (3.3)-(3.5) to the case with vanishing initial data and an inhomogeneous 
term that vanishes to all orders as t —► 0, by replacing W by W — W and F by F — LW in (3.3). 

Let us now introduce some notation: 

Definition 3.1. Let 

(3.12) \\W(t)\\ H r = \\W(t, Oilmen), 
and 

(3.13) 1|wm=£ \\b k W{t)\\ H s. 

s+fc<r 


Let N be the exterior unit normal to dD in the metric g a b, or equivalently, N a = N l dy a /dx l . Set 
(3.14) (W(t)) r = \\W N (t,-)\\ Hr{an) . 

where Wn = W ■ N is the normal component. 
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Theorem 3.1. Suppose that p = p(p) is a strictly increasing smooth function of p. Suppose also that 
x is a smooth solution of (2.7), such that (1.9) hold for 0 < t < T. Suppose that the inhomogeneous 
term F in (3.3) is smooth for 0 <t<T. Suppose also that the initial conditions (3.f) are smooth and 
satisfy the m th order compatibility conditions (3.8), for all m = 0,1,.... Then the linearized equations 
(3.3)-(3.5) have a smooth solution for 0 < t < T. 

Let 


(3.14) E r (t) = \\W(t)\\ r + \\W(t)\\ r + (W(t)) r + \\divW(t)\\ r , 


where W = D t W = D t W + (di vV)W. Then there is a constant C depending only on x, r and T such 
that for 0 < t < T we have 


(3.15) 


E r (t) < C(E r { 0)+ / \\F\\ r dr 


Theorem 1.1 follows from Theorem 3.1 since the norm (3.14) is equivalent to 

(3.16) E r (t) = \\W{t)\\ H r + \\W(t)\\ H r + (W(t)) H r + || divW(t)\\ H r, 

if F vanishes. In fact, by (3.6) one can express time derivatives in terms of space derivatives of the 
same order or less and using induction it follows that 

(3.17) E r < E r F C r (E r + ||J , | r _i) 

In this section we show the lowest order energy estimates for an equation of the form 

(3.18) W + CW = B(W,W) + F 

where W = D t W = k^D^kW), W = D(W, 

(3.19) 

CW a = -g ab d b (p' (divW + {d c e)W c )) = -g ab d b (h'div{pW)), e(p) = In p, p'{p) = h'{p)p 

and B is any bounded linear operator. The energy is: 

(3.20) E = {W,pW) + (W,p{C + I)W) 

= f g ab W a W b + g ab W a W b + p' (div (p W) / p) 2 pndyF f W^(-V N p)dS. 

Jn Jdn 

Now, for any symmetric operator B we have 

(3.21) —( W,pBW) = — / nW a pB_W a dy = 2(W,BW) + ( W,pB'W ), 

dt dt J n 

where W = D t (nW) and pB' is the time derivative of the operator pB considered as an operator 
from the vector fields to the one forms: 

pB'W a = g ab (D t ( P BW b ) - pBW b ), 
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(3.22) 


BW b = g bc BW c , 



Since ( W, pW) = (W,pGW), where G = I, it follows that 

(3.23) E = 2 (W, pW + p(C + I)W) + (W, p G'W) + (W, p (C' + G')W). 

where pG'W a = D t (p g a bK)W b and C/W a = D t CW a — CW a + eCW a , where e = p/p. Since D t (pn) = 0 
and D t (ndiv (pW)) = Kdiv (D t (pW)) we get 

(3.24) D t CW a = -d a D t ((pK)- 1 p , (p)Kdiv(pW)) 

= -d a (p"(p) e di v(pW) +p'(p)p~ 1 div (p W + pW)^j 
so 

(3.25) Q/W a = -d a (p"{p)e di v(pW)) - d a (p\p)p~ 1 div(p W)) - ed a (p'{p)p~ 1 div(p W)) 

Since p\ an = 0 and e = p/p it follows that 

(3.26) (■ U,pC'W) 

pp" div (p U) div (pW) + p'(div(pU) div(pW) + div(p 17) div(p W))^ p~ 1 Kdy 

+ [ ( -V N p)U N W N dS. 
Jon 

It therefore follows that E < C\[E(\[E + ||F||) and hence with E 0 = \[E we have 

(3.27) E 0 (t) < C 
4. Decomposition of the linearized equations into an operator on the 

DIVERGENCE FREE VECTOR FIELDS AND AN OPERATOR ON THE ORTHOGONAL COMPLEMENT. 

We will now make an orthogonal decomposition: 7i = L 2 = TLq © 7i\ into divergence free vector 
fields TLq and gradients of functions in Hq(Q); TL\. Let us therefore define the orthogonal projection P 
onto divergence free vector fields by 

(4.1) PU a = U a -g ab d bPu , A Pu = div U , Pu \ m = 0 

(Here A q = K~ 1 d a (ng ab dbq)■ ) P is the orthogonal projection in the inner product, see [LI], 

(4.2) (U,W)= f g ab U a W b ndy. 

Jn 

Note also that, with 

S 

(4.3) \mt)\\ r , s = ^\\D^W(t,-)\\ H r 

k =0 

denoting the Sobolev norms for fixed time with space and time differentiation of order r and s, we have 

(4.4) \\PW\\ r , s < C\\W\\ r , s , ||(L — P)W\\ r , s < C\\ divW|| r _ 1)S 
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since its just a matter of solving the Dirichlet problem and commuting through time derivatives, [LI]. 
For a function / that vanishes on the boundary define AfW = —P(y(W c d c f))', i.e. 

(4.5) A f W a = —g ab d b ((d c f)W c — q), where A((d c /)IT C - g) = 0, q\ dn = 0 
If U and IT are divergence free then 

(4.6) (U, A f W) = [ n a U a (—d c f)W c dS 

JdQ. 

where n is the unit conormal. If f\ tm = 0 then —d c f | SQ = (—Vj vf)n c . It follows that Af is a symmetric 
operator on divergence free vector fields, and in particular the normal operator 

(4.7) A = A h , 

where h is the enthalpy, is positive, i.e. (II 7 , AW) > 0, since we assumed the physical condition that 
—Vjy/i > Co > 0 on the boundary. The normal operator is order one, by (4.4) 

(4.8) \\AW\\ r , s < Cj|IF|| r+liS . 

The normal operator has certain delicate commutator properties with vector fields and positivity 
properties which were essential for the existence proof in [LI]. The main difficulty being that it is not 
elliptic acting on vector fields with non vanishing curl. In order to prove existence one had to replace it 
by a sequence of bounded operator which uniformly had the same commutator and positivity properties. 

We now make the decomposition 

(4.9) W = W 0 + W U W 0 = PW e Ho, W 1 = (I-P)WeH 1 
We want to decompose the linearized operator 

(4.10) LW = W + CW-B(W,W), divW\ m = 0 
where B is a bounded operator and 

(4.11) CW a = -g ab d b ((d c h)W c + p' divW), 

into an operator onto the divergence free part and an operator on the complement. 

The projection to highest order commutes with time differentiation: 

(4.12) PW 0 = W 0 , PW! = PB 2 (W 1 ,W 1 ), B 2 (W,W) a = -g ab (g bc W c + 2g bc W c ) 

where g ab = D t g a b, 9ab = D^g ab , and D t = D t -(divV). In fact, applying L> t 2 to g a bW b = d q q x gives 
9 abWi + 2g ab W b + g a bW b = d a qi. Here qi = D^qi vanishes on the boundary since (p does. The pro¬ 
jection of g ab d b qi therefore vanishes and (4.12) follows since D t preserves the divergence free condition. 
Furthermore with A given by (4.5)-(4.7) and C by (4.11) we have 

(4.13) PCW = AW, if div W | m = 0, 

since the projection of the highest order term, V [p 1 div W ), vanishes since div IT| = 0. 

We now want to project (4.10) onto the divergence free vector fields using (4.12)-(4.13). We get 

(4.14) PLW = W 0 + AW 0 + PB 2 (W 1 ,W 1 ) + PAW X - PB (II 7 , IT) 
where A is the normal operator (4.7). Similarly, applying (I — P) to (4.10) gives 

(4.15) (/ - P)LW = ITi - PB 2 {W 1 ,W 1 ) - V(p'divITi) - (I- P)V(W c d c h ) - (I- P)B{W ,IT), 
subject to the boundary condition div IT 1 1 dQ = 0. 
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Lemma 4.1. Let L be defined by 


(4.16) Lw 0 = W 0 + AW 0 , 

(4.17) LW 1 = Wi — PB 2 (W 1 ,W 1 ) - V (p' cliv W \) 

and let M be defined by 

(4.18) PMW = PB 2 (W 1 ,W 1 ) + PAW X - PB(W , IT), 

(4.19) (/ - P)MW = -(/ - P)V(W c d c h ) - (7 - P)B(W, W) 

We have 

(4.20) LW = LW + MW 


If Po = P, P\ = (/ — P) and Lij = PiLPj then L respectively M are essentially the diagonal 
respectively the off-diagonal part of L. It turns out that we can invert (4.16) on TLq, see section 5, and 
(4.17) on 7i\, see section 6. The interaction term M is lower order, but in subtle way, since it contains 
space derivatives dW. The estimates for (4.17) gives us control of an additional space derivative of 
W\ and that is all that is needed to estimate (4.18). (4.19) also contains a space derivative of Wo- 
However, in our estimates for (4.17) we can replace this space derivative by a time derivative and the 
estimates for (4.16) gives us control of an additional time derivative. The estimates for (4.16)-(4.17) 
will be summarized in section 7. 

5. Existence and estimates in the divergence free class. 

In [LI] we proved existence of solutions for 

(5.1) W 0 + AW 0 = F 0 , Wo| t=0 = Woo, W 0 \ t=0 = W w 

We have 

Proposition 5.1. Suppose that x,h are smooth, h\ gn = 0 and V/v/t < —Co < 0 for 0 < t < T. 
Then if initial data and the inhomogeneous term in (5.1) are smooth and divergence free it follows that 

(5.1) has a smooth solution for 0 <t<T. Furthermore, with a constant C depending only on the norm 
of x and h, T and the constant Co we have 

(5.2) E 0r (t) < c(^E Or (0) + J ||F 0 (T)||^dr), E 0r (t) = \\W 0 (t)\\ H r + \\W 0 (t)\\ Hr + (W 0 (t)) r 
where 

(5-3) \\W(t,r)\\ H r = \\W(t, ■ )||ij’-(n)i (W(t)) r = \\W N (t,-)\\ Hrm) 

and Wn = N a W a is the normal component. 

Proof. In case div V = 0 this was proven in [LI] and the proof there can be easily modified by multiplying 
or dividing by R = det (dx/dy). Let us now indicate what needs to be changed in [LI] in order to deal 
with the case divH^O. We can use the same set of tangential vector fields as in [LI], but they are no 
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longer divergence free so the Lie derivative with respect to these no longer preserves the divergence free 
condition. But one can easily modify the Lie derivative so it preserves the divergence free condition. 
The modified Lie derivative with respect to a vector field T applied to a vector field W is 

(5.4) C T W = C T W + (divT)W, 

where CjTV is the Lie derivative. It satisfies divT^IT = TdivIT, where for a function /, Tf = 
Tf + (di vT)f. One then has to make it so one always apply this modified Lie derivative to vector 
fields. However we use the usual Lie derivative, when applied to one forms since it commutes with 
covariant differentiation. In deriving the estimates for all components of a vector field in terms of the 
divergence, the curl and the tangential components we use jCT(dabW b ) = g a b^rW b + (jC.T 9 ab)W b , where 
-C-Tdab = C'Tdab— (divT) Let us now examine how the critical commutator with the normal operator, 

(4.5) , is changed from what it was in [LI]. With Af given by (4.5) and AfW a = g a bAfW b we have 

(5.5) C T AfW a = C T d a {(d c f)W c - q) = d a ((d c f)C T W + 3 c {Tf)W c -Tq + f{d c d^T)W c ) 

where Tf = Tf — (di vT)f. When we project again the last two terms vanish since they vanish on the 
boundary so the commutator relation in [LI] will be replaced by 

(5.6) PC T A f W = A f C T W + A ff W 

The issue of how to deal with the initial conditions in case divH A 0 was discussed in section 3. □ 

Now, the norms used in Proposition 5.1 are natural for the initial value problem. However, when 
solving the wave equation with Dirichlet boundary conditions it is more natural to first look on norms 
with many time derivatives. Because of the coupling between the two equations we must therefore also 
estimate more time derivatives of the divergence free part. From differentiating (5.1) we get: 

Lemma 5.2. Suppose that W is a smooth solution of (5.1) for 0 <t <T. Let E (i r be as in (5.2) and 

S 

(5.7) \\W(t)\\ r = \\W(t)\\ k , s where || W(t)\\r,s = J2 H ^W(t)\\ H r 

s+/c<r k =0 

Then for 0 < t < T we have 

(5.8) j|^o||r + ||Wo|| r <C7(£?o,r + ||i ; b||r- 1 ), T>\ 

Proof. The proof is just differentiation of (5.1) using that A is order one, (4.8); 

(5.9) HWollr—r <C(||Wo|| r + ||Fo|| r —0, r>l 
which proves (5.8) for r = 1 so we may assume that r > 2 in (5.8). We have 

(5.10) ||Wb|| r <||Wb|| r> o + ||^o||r-i,o + ||^o||r-2, r>2 

so together with (5.9) we get 

(5.11) 11 IFo 11 r < C(E 0 , r + HH'ollr-i + ||7'o ||i— 2 ) j r > 2 

Since also |ITo|| 1 < E 0 \ we can use induction in r to prove that 

(5.12) j|Wo||r < C(^o,r + ll^ollr- 2 ), r> 2 

Similarly 

(5.13) HWollr < ||Wo||r-,o + HWollr-— 1 , r>l 

so by (5.9) again 

(5.14) \\Wo\\r<C(Eo,r + 11^0^ + 11^011^!), T>\ 

(5.8) for r > 2 now follows from (5.12) and (5.14). □ 
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Theorem 5.3. With notation and assumptions as in Proposition 5.1 and Lemma 5.2 we have 


(5.15) E 0 r(t) < c(E Or (0) + ||F 0 (t)|Ut), E 0r (t) = \\W 0 (t)\\ r + ||Wo(t)|| r + {Wo(t)) r 

Proof. (5.15) follows from Proposition 5.1 since 

(5.16) ||F 0 (f)|| ? ._ 1 < ||F o (0)|| r _ 1 + / ||To(r)|| r ._ 1 dr , ||F 0 (0)|| r _ 1 < C(||lT(0)|| r + \\W(0)\\ r ) □ 

Jo 

6. Existence and estimates for the wave equation. 

We consider the Cauchy problem for the wave equation on a bounded domain with Dirichlet bound¬ 
ary conditions: 

(6.1) D 2 (e'if)-Aif = f, in [0,T] x fi, if\ gQ = °> 

(6-2) ^\t=o= Dt M t=o = ^ 

Here 

(6.3) A if = -^= d a ^ \/ det gg ah df/ifj , 

where g ab is the inverse of the metric g a b and det g = det{g a b} = k 2 , in our earlier notation. We assume 
that g ab is symmetric (since the metric is), and that g ab and e' are smooth satisfying: 

n 

(6.4) 0 < e'+ 1/e'< ci, ^ {\g ab \ + \g a b\) < nc\ 

a,6=1 

for some constants 0 < c\ < c' x < oo. 

Existence of solutions for a wave equation with Dirichlet boundary conditions and initial conditions 
satisfying some compatibility conditions is well known, see e.g. [H,Ev], In order for (6.1)-(6.2) to be 
solvable initial data must be compatible with the boundary condition. If we move the Laplacian in 
(6.1) over to the right hand side and differentiate (6.1) with respect to time we get 

(6.5) D^ 2 ' t f = b k ^,...,Df +1 ^d^...,dD^,d 2 ^...,d 2 D^)+D^f 
for some functions bk- We therefore define functions of the space variables only ipk 

(6.6) i> k+2 = b k (if o,..., ij>k+i,dij>o ,..., dip k ,d 2 ip 0 ,..., d 2 if k )\ f=0 + D'ff | (=Q 

where if o and if\ are as in (6.2). For this to be compatible with the boundary conditions we must have 

(6.7) V’fclan = 0) f° r k <m - 1 

(6.7) is called the m th order compatibility condition. Since if k are determined from the initial conditions 
ifo and ifi this gives some compatibility conditions on the initial conditions. We have: 

20 



Proposition 6.0. Suppose that g,e' are smooth satisfying ( 6 . 4 )■ Then if initial data (if 0 ,^ 1 ) and f 
are smooth and satisfy the m th order compatibility condition for all m, it follows that (6.1)-(6.2) has a 
smooth solution if. 

Proof. The result in [H] is stated with vanishing initial conditions but, if the compatibility conditions 
are satisfied, one can reduce to that case by subtracting off an approximate solution satisfying the 
equation to all orders as t —> 0. Let 

OO 

(6.8) = E X{t/£k)t k 'tik/k\, 

k =0 

where x is smooth y(s) = 1 for |s| < 1/2 and x( s ) = 0 for |s| > 1, and the sequence £& > 0 are 
chosen small enough so that the series converges in C m ([0,T], H m ) for any m. This is obtained if take 
(IIV’fc||j?' c + l) £ fc < 1/2. Then ip = if — if satisfies (6.1) with vanishing initial conditions and a right hand 
side / that vanishes to all orders as t — * 0: 

(6.9) AVA - = - AVA + + / = 7 

For this case existence of a smooth solution if to (6.9) follows from Theorem 24.1.1 in [H]. Since the 
theorem in [H] is more general let us just point out the main steps needed for our case. Existence 
follows from duality, using the Hahn-Banach extension theorem and the Riesz representation theorem. 
For this one has to show estimates for the adjoint operator in negative Sobolev spaces. Suppose that if 
satisfy (6.1) and let ifw = (I — A)~ N if, where N > 0 and A is the Dirichlet Laplacian, i.e. inductively, 
we define ifk to be the solutions of (I — A)ifk+i = ifk, with boundary conditions ifk+i\ i)n = 0. Then 
ipN satisfies (6.1) with / replaced by + (I — A)~ N [D^,(I — A) A ]ipN, where /jv = (/— A)~ N f. The 
norm of this is bounded by ||/at|| + IIVAI + lAt^jvll- Using the energy estimate in Lemma 6.2 then 
gives us an estimate -)|| + \\X7ipiv(t, -)|| < C f ( ' || /at || dr □ 

Lemma 6.2. Suppose that g ab ,e\f and if are smooth and satisfy (6.1)-(6-4) for 0 < t < T. Let 
V“ = g ab db and, for r > 1, 


( 6 . 10 ) 


1 -1 


m = E 


s=0 


e\D t s+ 1 if) 2 +\D s t ViP \ 2 + if‘ 


) ndyj 


1/2 


Then 

(6.11) — <CA+||/||o, r -i), where U\\ r ,s = E W D t d y < t > II 

fc<s, |a|<r 


Proof. We will prove that dE 2 /dt is bounded by E times the right hand side of (6.11) and (6.11) 
follows from this since dE/dt = (dE 2 /dt)/(2E). Since D t n = K.divU, we have with D t = D t +divU 
and Dt = Dt — div V: 

( 6 . 12 ) E ^(e'(A + V)(A + 2 A+ 5 a 6 (Av a A(A + 1 V b A)«^ 

+ \ f ((Ae')(A + V) 2 + (■ D t g ab ){D S tV a if) {D?V b if) + K~\D t K)if 2 + 2 ifD t if) ndy. 
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Here the terms on the second row are bounded by a constant times E 2 . Applying D( +1 to d a ^ = g a b'^ b ' t l ) 
gives d a D t s+ V = E*=o (E) (^t +1 ~ % 9ab)DiV b (f> so 


(6.13) 


A S+1 V a V = 9 ab d a D s t + ^ - J2 ( a+ i 1 )9 ab (D a t +1 ~ i 9 bc)Di^- 


i =0 


Up to terms bounded by a constant times E 2 , (6.12) is therefore equal to 

(6.14) / (e'(A + V)(A + A) + (A'V» (a o D-+V)) «<& 

s<r— 1^^ 


E f ((A + V)(e'A + V - wd*/ 

s<r 


where we have integrated by parts. If we apply £>f to D 2 (e'if) — k 1 d a ( y nV a if) = f we obtain 


s+l 


(6.15) A(A(eV) - /^(/eV^)) = e'A + V - « _1 3a(4 s V» + £ ( s f)(A +2 “V)(AV0- 

i=0 

Since the L 2 norm of the last term is bounded by CE plus the L 2 norm of if the lemma follows. □ 

One can get additional space regularity from taking time derivatives of the equation (6.1) and 
solving the Dirichlet problem for the Laplacian. 

Lemma 6.3. Suppose that g ab and e! are smooth and satisfy ( 6 . 4 ) and that f is smooth, for 0 < t < T. 
Suppose also that ip is a smooth solution of (6.1), for 0 <t <T. Let ||^||' S}T . he as in Theorem 6.1 and 
HWr = E s+fc <r \\ D (M\k- Then 


(6-16) HV’llr < C(\\lf\\ 0 , r + HV’llgr-l + ll/ll,— 2 ) - 

Proof. Since A if = Z) 2 ( e V) — / and if\ gn = 0 it follows that 

s+2 

(6.17) A D\$ = E ( a t 2 )(A 8 + 2 _i e')A^ - D?f ~ [4 S A - A D?}if D$if | fln = 0 


i =0 


so by the standard elliptic estimates 


s -|-2 


s —1 


(6.18) \\D^i/j\\ H k+2 <c(E« T>lip\\H k + II T>(f\\ H k + E \\D l t 'ip\\ H k+2) 

i =0 i —0 

Here the last term is lower order and is absent if s = 0 so using induction in s we get 


(6.19) 


EIIAVI \hw < cJ2 (IIA + Vll^ + II D\f\\ Hk ) 


»=0 


i =0 


or with ||V>||r,» = EL 0 W D tV\\H' 

(6.20) ||V’|| s +2,r- s -2 < C(W\., r - a + \\f\\s,r- s -2), 0 <s<r-2 

Since ||V’||r = El=o ||V’l|r-a,s it therefore inductively follows that 

(6.21) ll^llr < C'dl^llo.r + II^H l,r —1 + ||/||,- 2 ) □ 


Summing up; 
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Proposition 6.4. There are constants C r such that the solution of (6.1) satisfy 


( 6 . 22 ) 


IIVAOIIr < C r (\\ip(0, -)||r- + 


-2 dr), r >2 


Proof. It follows from Lemma 6.2 that E(t) < C r (E( 0) + ||/||o,?—l dr). Using (6.13) we see that the 

energy energy E in (6.10) is equivalent to ||V ; ||o,r + ||V’l|i,r-i- Furthermore || f(t, -)||r -2 < 11/(0, -) 11 2 + 
/o || f(t, •) ||r —2 dr and since also ||/|| r _2 < C'||'0|| r the proposition follows from Lemma 6.3. □ 

As pointed out in section 4 we actually want to solve the equation: 

(6.23) m - PB 2 (W 1 ,W 1 ) - V(p'divWi) = F u divlUi| an = 0 

(6.24) W\ = Vqi, gi| aQ = 0, 

where (/ — P)F\ = F\, which is equivalent to 

(6.25) A 2 / - A(//r/>) = divA, /|an = 0 

(6.26) Wi = Vq 1 , Aqi = <j>, q x | an = 0. 

In fact, for W\ of the form (6.24) the left hand side of (6.23) is (I — P)W\ — (I — P) V (// div W \), and 
(/ — P)H = 0 is equivalent to di vH = 0. Assuming that the compatibility conditions are satisfied we 
can solve (6.25)-(6.26) and this then also gives us a solution of (6.23)-(6.24). The initial conditions for 
(6.25) are 

(6.27) f >\ t=0 = 4> 0 , D t (j) | t=Q = A 
Lemma 6.5. Suppose that W\ satisfies (6.23) and set 

(6.28) P(i) = (^E f {ID^W ^ 2 +p'\div(D s t W 1 )\ 2 PlW^Kdy) 

Z s=0 dn 


1/2 


Then 

(6.29) 


— <C{E+\\F 1 \\^ r ). 


Proof. Let W lk = DfW u 
dE 2 _ _ 

s =0 Jn 


(6.30) 


dt 


E/ {W ls -Wu+p'div W ls div W ls ) ndy 


+ \ f (( D t g ab )W? s Wl + (D t p')(dwW u ) 2 + K - 1 (D t K)\W 1 \ 2 + 2W 1 -W 1 )Kdy. 
1 Jn 

Integrating by parts, using that div!Ui s | an = 0, it therefore follows that 

dE 2 


(6.31) 


dt 


< 


E L W u • (W l8 - V(j/divIU l6 .)) Kdy + CE 2 
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s=0 



Using (6.23) this proves (6.29) for r = 0. To prove it for r > 1 we have to commute time derivatives 
through (6.23), which can be written 

(6.32) g ab W b -d a (p'divW 1 ) =g ab Bl B 0 (Wi,Wi, F\) = F 1 + PB 2 (W 1 ,W 1 ) 

With q = p' / k we have 

(6.33) D t s (g _1 a o (p , divWi)) = D)d a (ndwWi) + D s t ((d a lng)/«divWi) 

S S — 1 

= d a (KdwW ls ) ( s k )(d a D%~ k lnq)ndivW lk = q~ 1 d a (p / dwW ls ) + ^ (£) (d a D)~ k In q)nddvW lk 

k= 0 k= 0 

Hence 

S — 1 

(6.34) ^(g-^^p'divWi)) = d a (p'dwW ls ) + WaD s t ~ k lnq)p'dwW lk 

k= 0 

Multiplying (6.32) by g _1 , applying D) and dividing by g _1 therefore gives 

(6.35) g ab W b ls - d a (p'divW ls ) = g ab B b s (W 1 , ..., Wu.divW^, ....div^r, F u ..., DfF^ 
where B s is a bounded operator of its arguments: 

(6-36) ]|5 fc || <C^(||W lfc || + ||W lfc || + \\D k Fi || + ||divW lfc ||) 

k= 0 

This together with (6.31) proves (6.29) also for r > 1. □ 

Theorem 6.6. Suppose that the initial conditions (6.27) and the inhomogeneous term in (6.25) are 
smooth and satisfy the compatibility conditions for all orders. Then (6.25)-(6.27) has a smooth <f>. 
Furthermore with W\ given by (6.26) we have 

(6.37) E lr (t) < C r E lr (0) + C r [ ||T 1 (r)|| r _i dr, E lr (t) = ||Wi(f)|| r+ i 

Jo 

for r > 1, and for r = 0 the same inequality holds with ||Fi(r)||,._i replaced by ||Fi(r)|| r . 

Proof. First we assume that r > 2. By the second part of (4.4) we see that with <f and W\ as in (6.37) 

(6.38) ||0||r < \\9WiWr + \\WxWr < C\\(j>\\ r 
where d stands for space derivatives only. Furthermore by (6.23) 

(6.39) HfUilUi < WdW^l + ||Will,. + ||Wi|| r _i < C l (||(91U 1 || r + H^H,. + ||F 1 || r _ 1 ) 

Since also 

(6.40) ||Ti(t)|| r _i < ||Fi(0)|| r _i + f UTiHr-i dr < C , ||H'i(0)|| r _)_i + f ||Fi|| r _i dr 

Jo Jo 

(6.37) for r > 2 follows from Proposition 6.4. For r = 1, (6.37) follows from Lemma 6.5 and the fact 
that by the second part of (4.4) ||<9LhfUi|| < C(\\D t divVFi|| + ||divVFi||) and by (6.23) ||<9 2 fUi|| < 
C'HadivWiH < C'OlWiH + \\D t dWi || + ||W||i) + ||Ti||. □ 
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7. The proof of the theorem. 


We are now in position to prove Theorem 3.1. We want to show that 

(7.1) LW = F, divW| aQ = 0, 

with initial conditions 

(7.2) w | (=0 = w„, ir | t=0 = w„ 

has a smooth solution W if Wq , W\ and F are smooth and satisfy the compatibility conditions in 
section 3 to all orders. If these compatibility conditions hold then we can find an approximate solution 
W satisfying the initial conditions and the equation to all orders as t —> 0. Subtracting off this 
approximate solution reduces it to finding a smooth solution to (7.1) when Wo = W\ = 0 and F 
vanishes to all orders as t —» 0. 

With L and M are as in Lemma 4.1 we have reduced to finding a smooth solution of 

(7.3) LW = F, div W | gn = 0, 

where 

(7.4) F = F - MW. 

That the operator L is invertible follows from using the decomposition W = Wo + W\, in Lemma 4.1 
and applying Theorem 5.3 to the divergence free part Wo and Theorem 6.6 to W\. This in fact gives 
estimates for the solution of (7.3) that can be used to show existence for (7.3) with F given by (7.4) by 
iteration, and we hence obtain a solution to (7.1). Let us introduce the norms: 

(7.5) \\\W\\\ ril = ||Wb||r + ||Wo|| r + {Wo)r + HWdlUr, W 0 = PW, W 1 = {I-P)W 
and 

(7.6) |||T|||. r , 2 = ||Fo|| r + HAllr- 1 , Fo = PF, F 1 = (I — P)F. 

It now follows from Lemma 4.1, Theorem 5.3 and Theorem 6.6: 

Theorem 7.1. Suppose that (x,h) is a smooth solution to (2.7), (for 0 <t<T), such that h\ gn =0 
and V/v^lan < — Co <0. Suppose also that F, Wo and W\ are smooth and such that there is a smooth 
function W satisfying the initial conditions (7.2), the boundary condition div!T| an = 0 and LW = F to 
all orders as t— >0, i.e. [LW— F) 1 1=0 = 0, for k> 0 Then (7.2)-(7.3) has a smooth solution W. 
Furthermore, there are constants C r such that for any smooth solution of (7.3) we have 

(7.7) \\\W(t)\\\r, 1 <a(|||^(0)|||r.l+^ \\\F\\\r, 2 dTy V > 1. 

Moreover, 

\\\MW\\\ r , 2 < C r \\\W\\\ rA . 
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(7.8) 



We remark that the compatibility conditions in the theorem are in particular true if Wo = W\ = 0 
and D+F I „ = 0, for k > 0, since then we can take W = 0. Therefore if D+F\_, „ = 0, for k > 0 and 

6 I t=u L 1£=U 

we set W° = 0, it follows that we can inductively solve, for k > 0: 


(7.9) LW k+1 = F — MW k , divW k+1 \ m = 0, W k+1 \ t=Q = W k+1 \ t=0 = 0, 


We claim that W k converges to a solution of (7.3)-(7.4) and hence to (7.1), in case D k F\ f (| = 0, 
for k > 0 and Wo = W\ = 0. Since we have already reduced solving (7.1)-(7.2) to this case this 
would prove the existence part of Theorem 3.1. That W k converges to a smooth solution of (7.3)- 
(7.4) follows from using the estimate in Theorem 7.1. In fact, L(W X — W°) = F, and for k > 1 
L(W k+1 — W k ) = —M(W k — W k ~ 1 ). It therefore follows from Theorem 7.1 that 


r t n 

(7.10) M n <C r (|||T||| r)2 + Mjv) dr, where M N = ^ |||W fe+1 - W 

*^° k =0 


r, 1 


It now follows from a standard Gronwall type of argument that 


(7.11) 


Mw(t) < C r e 


c r t 


IFI 


\r, 2 


dr 


for any N. It follows from this that W k converges to a smooth solution W of (7.3)-(7.4) for 0 < t < T 
and therefore we have proven existence of smooth solutions for (7.1)-(7.2). 

Haven proven existence of a smooth solution to (7.1)-(7.2) we now also need to prove the estimate 
in Theorem 3.1. Applying the estimate in Theorem 7.1 gives 


(7.12) \\mm\r,i<c r (\\\wm\\r,i+£ limil^ + lllHlI^dr). 

By the same Gronwall type of argument as above we get 


(7.13) lllW^HIr,! <ae C ^(|||lT(0)||| r , 1 +^ |||F||| r , 2 dr), r > 1. 

It therefore only remains to observe that the norms in Theorem 3.1 are equivalent to those here. It 
follows from the continuity of the projection (4.4) that 

(7.14) ||W|| r + \\W\\ r < ||Wo||r- + HWillr- + ||Wo|| r + ||Wi|| r < a(||^||r + ||W||r) 

Furthermore, by the second part of (4.4) and Sobolev’s lemma 

(7.15) || divVF|| r < ||VFi|| r+ i < C , r ||divW'|| r , and {W x ) r < C r ||VFi || r+ i 


This concludes the proof of Theorem 3.1. 

Acknowledgments. I would like to thank Demetrios Christodoulou, David Ebin and Kate Okikiolu 
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